Formation of Non-Abelian Monopoles Connected by Strings 
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We study the formation of monopoles and strings in a model where SU(3) is spontaneously broken 
to (7(2) = [SU(2) x f7(l)]/Za, and then to 1/(1). The first symmetry breaking generates monopoles 
with both SU(2) and t/(l) charges since the vacuum manifold is CP 2 . To study the formation of 
these monopoles, we explicitly describe an algorithm to detect topologically non-trivial mappings on 
CP 2 . The second symmetry breaking creates 1i strings linking either monopole-monopole pairs or 
monopole-antimonopole pairs. When the strings pull the monopoles together they may create stable 
monopoles of charge 2 or else annihilate. We determine the length distribution of strings and the 
fraction of monopoles that will survive after the second symmetry breaking. Possible implications 
for topological defects produced from the spontaneous breaking of even larger symmetry groups, as 
in Grand Unified models, are discussed. 



Topological defects are formed in a vast array of labo- 
ratory systems and may also have formed during a cos- 
mological phase transition The statistical proper- 
ties at formation of the simplest of defects have been 
studied quite extensively in the context of cosmology @ 
and more recently in a variety of different condensed- 
matter systems. Experiments have been performed to 
observe the spontaneous formation of defects in nematic 





liquid crystals 0. 
superconductors 




m superfluid 3 He % 0] and in 
In most particle physics applica- 
tions, the vacuum manifold can be quite complex, and 
hybrid topological defects may be formed. These may 
consist of monopoles connected by strin gs o r walls that 
are bounded by strings (see for example |lG|). 

In this paper we study the formation of non-Abelian 
monopoles that subsequently get connected by strings 
due to a second non-Abelian symmetry breaking. More 
specifically, we study monopoles formed in the symmetry 
breaking 



SU(3) -> 17(2) ee [SU(2) x C7(l)]/Z 2 . 



(1) 



The fundamental monopoles carry both SU(2) and £7(1) 
charge and may be labeled by a pair of charges, (1, ±1), 
where the first entry (with no sign) is the SU(2) charge, 
and the second entry is the 1/(1) charge. After the 
monopoles are formed, we consider the further symmetry 
breaking 



SU(2) ->Za. 



(2) 



Now all the monopoles will get connected by strings. 
However, the SU (2) charge is a Z 2 charge, and so there 
are two types of monopole states connected by strings 
(Fig. [TJ. The first of these is a monopole-antimonopole 
bound state i.e. a bound state of (1, +1) and (1, —1). The 
confining strings will then eventually bring the monopole 
and antimonopole together and lead to their annihila- 
tion. The second possibility is that the string confines a 
monopole to a monopole i.e. two (1,+1) or two (1,-1) 
objects. In this case, the confining string will bring to- 
gether the two monopoles to form a charge 2 object, 



FIG. 1: Two types of confined monopoles in the SU (3) model. 
The picture on the left represents a monopole and an anti- 
monopole connected by a string. The picture on the right 
shows two monopoles with the same C/(l) charge connected 
by a string. 



(0, ±2), that carries no net SU (2) charge but carries twice 
the basic £7(1) charge. One of our aims is to determine 
the relative number densities of the two types of objects 
subsequent to the second symmetry breaking stage. 

In the context of Grand Unification Theories (GUTs), 
fundamental magnetic monopoles also carry non-Abelian 
charges. For example, in the minimal GUT model 
with SU(5) symmetry, the fundamental monopoles carry 
51/(3) color, SU{2) weak, and U(l) hypercharge quan- 
tum numbers. The formation of magnetic monopoles in 
the grand unified context occurs due to the non-trivial 
topology of a very large vacuum manifold and our toy 
SU(3) model may be expected to capture some of the 
complications. 

One motivation for considering the formation of strings 
that connect non-Abelian monopoles is that the physics 
of confinement is not fully understood, and it is possible 
that non-Abelian magnetic fields also get confined due to 
quantum or plasma effects [IJ (HI ■ A second related mo- 
tivation comes from the Langacker-Pi proposal to solve 
the cosmic monopole over- abundance problem flU ] . The 
scenario assumes that electromagnetic gauge symmetry 
is spontaneously broken for a period in the early universe. 
As a result, magnetic monopoles carrying electromag- 
netic flux will get confined by strings and annihilate ef- 
fectively. Later the electromagnetic symmetry is restored 
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to be consistent with present observations. The breaking 
of SU (2) in our toy model performs a similar function for 
this non-Abelian model as does the Langacker-Pi mech- 
anism for the Abelian case, although it does not involve 
symmetry restoration at low energy. Monopoles again 
get connected by strings but here they can either annihi- 
late or form charge 2 states. The corresponding scenario 
in GUTs is more complicated since the monopoles get 
connected by several different kinds of strings [ll], G3 , as 
we discuss in Sec. [TV] 

We start in Sec. U by describing the field theoretic 
model under consideration, focussing on the topological 
aspects. In Sec. [II] we describe our numerical implemen- 
tation to study defect formation in the model and the 
results in Sec. IIIII We conclude in Sec. IIVI by discussing 
defect formation in an SU(5) GUT model. 



I. MODEL 

Our model contains an SU(3) adjoint field, <!>, whose 
vacuum expectation value (VEV) implements the sym- 
metry breaking in Eq. JT]). Two more SU(3) adjoint 
fields, *i and * 2 , acquire VEVs to break the SU(2) sub- 
group of U(2) to Z2 as in Eq. ©. The Lagrangian for 
the model is 

L = itr[(^$) 2 ] + i^tr[(^*0 2 ] 



i=i 



■-tr(VJH-^»*i^2), 



(3) 



where = 9^$ — ig[X^^ $], X^ LV is the field strength 
for the SU(3) gauge field V M , and the potential, V, is 
assumed to have a form that is suitable to give the fields 
the desired VEVs. 

The first stage of symmetry breaking is achieved by 
the VEV 



$ = $(°) = j-jT 8 




(4) 



where rj is the energy scale at which the first symmetry 
breaking occurs and will be set to unity since its value 
has no effect on the topological structures we are consid- 
ering. (We could also take $ = g&°^g^ for any global 
g G SU(3).) The vacuum manifold at this stage is 



SU{3)/U{2) = CP 2 



(5) 



Points on CP 2 are labeled by three complex numbers 
(zx, Z2, Z3), identified under a (complex) rescaling 

Z T = (z 1 ,z 2 ,z 3 ) = n(z 1 ,z 2 , z 3 ) , K € C, k ^ 0. (6) 

It will be convenient for us to label the points, following 
by a point on an octant of a two-sphere given by 9 
and <f>, and two phases, a and (3: 

Z T = (sin cos e M ,sin0sin0 e l/3 ,cos0), (7) 



with < 6, < tt/2 and < a, < 2tt. 

The relation between the field <& and a point on CP 2 

is 



1 / Ztf\ 
$ = -= 1- 3^— 



(8) 



The second homotopy group of CP 2 is known to be the 
set of integers Z. A topologically non-trivial configura- 
tion can be constructed explicitly by taking = 0. The 
points on the = sub-manifold are 



Z T = (smO e M ,0,cos< 



(9) 



and these describe a CP 1 subspacc of CP 2 . The points 
on a two-sphere in physical space, labeled by (0,0), can 
be mapped onto this CP 1 using 

= 0/2, = 0, a = <t>, (3 = 0. (10) 

Equivalcntly, 



1 



3cos0-l -3sin0e^' 

= fm I 2 

2 v3 V_3sin0e-^ -3cos0-l, 



This map represents a simple example of a monopolc. 

An expression for the topological charge of a monopolc 
can be derived by first constructing the 1-form "gauge 
potential" 



A = 



1 zUz - dz^z 



21 



Z^Z 



(12) 



Note that under the "gauge transformation" Z — > Ze lX , 
which is a special case of ©, A transforms as A — > A+dX. 
The corresponding field strength 2-form is 



, „ 1 fdZ^AdZ 



dZ^Z KZUZ\ 

(zTzy ) 



(13) 



Since this 2-form is exact, its integral over a closed two- 
surface is a topological invariant — and moreover is zero 
unless the surface contains in its interior a point or points 
where Z = (so that A is undefined). So the expres- 
sion for the topological charge in a volume V with closed 
boundary dV is 

Q=±- ( F=±- [ d 2 S*e^ k F jk . (14) 
2tt J dv 4tt J 9v 

There is another way to obtain the expression for the 
topological charge. We start with the expression known 
for the 't-Hooft-Polyakov monopolc in SU(2) and extend 
it to SU{3): 

Q = ±- f d 2 S l f abc e^ k n a d jn b d k n c , (15) 

Joy 



where 



Z^T a Z 
ZiZ ' 



(16) 
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with a,b,c = 1, . . . , 8. Here the T a are the generators 
of SU(S), normalized by tr(T a T b ) = 2S ab , the f abc are 
structure constants defined by [T a ,T b ] = 2if a b c T c , and 
the integration is over the two sphere at infinity. Also 
note that the vector n a satisfies n a n a = 4/3. In Ap- 
pendix [A] we show that the two forms for the topolgical 
charge are equivalent. 

It is simple to check that Q = 1 for the monopole 
configuration in Eq. (fT0|) and Eq. (fTTj) . The formula in 
Eq. p4p will be useful to locate monopoles in our numer- 
ical work described in Sec. [TX| 

The second stage of symmetry breaking is more in- 
volved. The fields now also acquire VEVs, which 
are required to lie in the unbroken SU(2) subgroup, and 
hence commute with Their magnitudes tr^ 2 ) are 
fixed by the potential, and they are also required to be 
mutually orthogonal in the sense that tr^i*!^) = 0. 
Given a value of $ at some spatial point P, we need 
to identify this unbroken subgroup. The standard proce- 
dure is to work out commutators of <E> with SU (3) gen- 
erators and to find linear combinations of the generators 
that commute. In practice, it is easier to first rotate 
$, say by an SU(3) rotation R, to the reference direc- 
tion, <J>(°). We discuss how to choose R below. Then 
the generators of the unbroken SU{2) sit in the 2x2 
upper left corner while the generator T 8 of the unbroken 
U(l) is in the direction of itself. With respect to 
the VEVs of *i and * 2 can be written in terms of 



two orthonormal 3-vectors, a and b, as $ 
*' 0) = b T where 



(0) 




a • T and 



(17) 



and Cj are the Pauli spin matrices. Once \I/ 1 ' ) and ^ 2 °^ 
are constructed, we can rotate all the fields back to the 
original point using RK 

The VEVs of *i and * 2 break SU{2) down to Z 2 , 
which is the center of SU(2), {1,— 1 2 }, i.e. the identity 
element of SU(3) and — 1 2 = diag(— 1, — 1, 1). A string 
passes through a spatial contour if and ^ 2 are such 
that, on going around the contour, these fields are trans- 
formed by the element — 1 2 and not by the identity el- 
ement. The strings are of the Z 2 variety and there is 
no distinction between a string and an anti-string. Also, 
there is no known integral formula that can be used to 
evaluate the winding around the contour. 



II. NUMERICAL IMPLEMENTATION 

To simulate the formation of the monopole-string net- 
work, a 3-dimensional cubic lattice is chosen. Each cu- 
bic cell is further divided into 24 tetrahedral sub-cells, 
obtained by connecting the center of the cube to the 8 
corners and the centers of the 6 faces (see Fig. [2]). 
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FIG. 2: Each cell of the cubic lattice is sub-divided into 24 
tetrahedra. Only one cubic cell and representative tetrahe- 
dron are shown. 



The next step is to assign random points of CP 2 at 
each point on the lattice, including the centers of the 
cubic cells and their faces. Now, the unique SU(3)- 
invariant metric on CP 2 is the Fubini-Study metric 



ds 2 



dzUz dz^zzUz 



z^z {z^zy 2 ' 

or, in terms of the parameter choice of l|7|. 
ds 2 = d0 2 + sin 2 6dd> 2 



(18) 



sin (9 cos 0(1 — sin #cos <f)da 



—2 sin 
,2 



' cos 2 </> sin </> da d(3 



sin- V sin 2 0(1 - sin 2 9 sin 2 0)d/3 2 . (19) 

Hence the S'[/(3)-invariant measure on CP 2 is 

y/g d9 dip da d(3 = sin 3 9 cos 9 sin <j> cos <f> d9 dip da d(3. 

4 - (20) 

Thus the assignment is done by drawing < sin 9 < 1, 
< sin 2 4> < 1, < a < 2-k and < /3 < 2ir from uniform 
distributions, and then constructing Z as in Eq. ([7]). The 
four vertices of a spatial tetrahedron then get mapped 
on to a tetrahedron in CP 2 which we will denote by 
{Z 1 ,Z 2 , Z 3 ,Z 4 ). To find out if this tetrahedron in CP 2 
is topologically non-trivial (i.e. incontractable) we use a 
discrete version of the charge formula in Eq. (fT4]) 



(21) 



{ijk} 



where the sum is over the four triangular faces of the 
tetrahedron (with positive orientation), and for each face, 



Y {ijk\ 



MgiZlZkZtZjZ^Zi), 



(22) 



where we require <%{ijk} to lie within the range [— 7T, it]. 
We can explicitly check that small changes in the Zi do 
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not affect Q, thus showing that even the discrete formula 
is topological. 

One can also check that Eq. (|2Tj) agrees with Eq. (fT4]) . 
The charge Q is the integral of F/2ir over a large sphere, 
which can be broken up into the sum of the four separate 
contributions from the individual faces of the tetrahe- 
dron. Each of these can be expressed as the integral of the 
1-form A/2tt around the perimeter. In discrctized form, 
the integral of A along the 1-2 link becomes axg(z\Z\) 
(see Eq. (|12p) and so the magnetic flux through the tri- 
angular plaquette {123}, is found by summing the con- 
tributions from the three edges, 



dx-A = arg(Z 2 t Zi) + arg(Z 3 t Z 2 ) + arg(Z 1 t Z 3 ) 

+2im, (23) 



so along a geodesic in CP 2 , i.e. RjiZi 



In 



fact the left-hand side is equal Zj times the phase 

factor that makes the scalar product with Z\ real 
(see Appendix IB"j) . In other words, we find 



RjiZi 



z zfa 

J \z}z t \ 



(25) 



The geodesic condition will be achieved if Rji can 
be written as 



R 



ji 



exp (iMs), 



(26) 



where M is a suitably chosen normalized combi- 
nation of the generators T a and s is the geodesic 
distance between Zi and Zj, given by 



where n is an integer and the extra term, 27rn, in Eq. (|23|) 
is included because each of the phases is ambiguous up 
to ±27r. This can also be seen as a gauge ambiguity: 
a gauge transformation may change the value of n. It 
has a geometric interpretation as well. For the special 
case of triangles on a CP 1 subspace of CP 2 (isometric 
to a sphere of radius 1/2), we have shown that the flux 
through a triangle, found using Eq. (fT4"|) . is equal to twice 
the area of the triangle. Thus the ambiguity in the flux 
in Eq. (|23p is equivalent to the ambiguity in choosing 
between the two complementary spherical triangles with 
this boundary. We choose the one with the smaller area, 
so that 



<j> dx ■ A 



*{123}- 



(24) 



Thus Eq. (|21j) is the discrctized version of Eq. (|T4j) . 

We conjecture that for a general triangle in CP 2 , not 
lying on a CP 1 subspace, the flux through it may still be 
equal to twice the area of the minimal surface with that 
boundary. Choosing the minimal area may be seen as a 
generalization to areas of the "geodesic rule" for lengths 
[2| . The rule in general is to choose the minimal value of 
the integral in Eq. (|23p . 

Next we turn to the formation of strings that connect 
the monopoles. For this we need to consider a triangular 
face of a tetrahedron and determine if a string passes 
through it. 

Each vertex of a triangular plaquette has already been 
assigned a point on CP 2 , equivalently a VEV of $. It is 
convenient to label the subgroup that leaves <&j invariant 
as SU(2)i x U(l) t /Z 2 . Now we also assign VEVs of *i 
and ^2, making sure that these lie in the unbroken SU (2) 
sector of SU(3) at Zj, namely SU(2)i, and that they are 
orthogonal: tr^i*!^) = 0. The precise scheme is as 
follows. 

• The scheme is based on the construction, for each 
pair of points on CP 2 , say Zi and Zj, of an SU(3) 
transformation, Rji, that transforms Zi to some 
representative of the point Zj and moreover does 



(ZjZj){zjZi) 
\\ {z\Zi){z\Zj] 



(27) 



A more explicit construction of Rji is described in 
Appendix IBl 

Similarly, for each Zj, we define an SU(3) transfor- 
mation Rio such that Zj = R^Zq, where Zq is the 
reference point (0,0,1). (With our choice of rep- 
resentative in ((?]), no phase factor is needed here.) 
The matrix R described in the previous section, 
above Eq. (fTT)) . will be one of the R' iQ . 

• To each vertex of the triangular face is associated 
a point on CP 2 (say Zj) and two uniformly dis- 
tributed orthonormal 3-vectors, aj and bj where i 
labels the vertex of the triangle (see Fig. [3]). If we 
wish, we can construct $j from Zj using Eq. I|8]). 
The two remaining fields ^x,2 may be found from 
a and b. We first define 



A i0 = a-T, P i0 =bT, 



(28) 



which are SU (3) matrices lying in the SU (2)o sub- 
group, with generators T given by Eq. (jTTJ) . Then 
the fields are given by \I>i = T)iA and ty 2 = V2B, 
where 771,2 are the magnitudes of these fields, and 
the normalized SU(3) matrices A and B may be 
found by using the transformation Pjo: 



Ai — RioAioR\ , Bi — PjoPjoPj, 



(29) 



Note that by construction Ai and Pj belong to 
SU(2)i and hence commute with $j. 

• Now we want to compare the symmetry-breaking 
fields at neighboring vertices. To do this we trans- 
port them using the geodesic transformations Rji. 
Transforming Ai and Pj by parallel transport along 
a geodesic from Zj to Zj, we obtain 



Aji — RjiAiRji 



Bji - /.',,/',/>',,. 



(30) 
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Next we compare these transported matrices with 
the corresponding matrices Aj , Bj defined at the 
vertex Zj. We seek a transformation Sji G SU(2)j 
such that 



Aj — SjiAjiSjt , 



Bj 



SjiBjiSji- 



(31) 



In Appendix [Cl we describe our construction of Sji 
in detail. 

The net rotation of the pair Ai , Bi as we circum- 
navigate the triangular face from Zi to Zj to Zj. 
and back to Zi is 



S 



{ijk} 



= .S',7-/l',/, .S'/, j I'i jS,, R 



(32) 



Note that since this combined transformation 
leaves invariant all the fields Ai, Bi, it must be- 
long to the unbroken U(l)i. 

• To determine whether or not a string passes 
through the {ijk} face, we have to compare Sujk} 
with the transformation RikRkjRji without the 
intervening S factors. Since this transformation 
leaves invariant, it belongs to SU{2)iX.U{l)i/1i2- 
Moreover, in view of Eq. (|22[) . we know that 



RikRkjRjiZi 



Consequently, wc know that the t7(l) 
product must be 



exp(- 



-,10L 



{ijk} 



(33) 

factor in this 



(34) 



Now let us return to S^jk}- Since for example the 
transformation Sji G SU (2)j leaves Zj unaltered, 
it is clear that, regardless of the choice of the S 
factors, the effect of S^jk} on Zi must be exactly 
the same as that of the product in Eq. (f3"3")) . Con- 
sequently, the combination 

W{ ijk } = S{ ijk y exp(±ia{ ijk }V3T?) (35) 

must leave Zi invariant, and also not contribute a 
phase to Zi, and hence it belongs to SU(2)i. But 
we know that W^k} also belongs to U(l)i, since it 
consists of two factors each of which is an element 
of U(l)i. So Wujfy must in fact be one of the two 
central elements that are common to both SU(2)i 



and U(l)i. If W, 



{ijk} 



1, the winding is trivial 



and there is no string through the triangular face. 
If, however, W^jk} = — l2j then there is a string 
through the triangular plaquette. 

It can be shown (see Appendix [D]) that if the 
monopole charge (|2"Tj) within the tetrahedron is 
non-zero, then there must be an odd number of 
faces with strings passing through, while if it is 
zero there must be an even number. This follows 
from the fact that each edge, say (ij) appears, with 
opposite orientation in two faces, and the relevant 
factors in say Sujfy and Syn} are inverses of each 



other: (SjiRji)^ 




b,) 



FIG. 3: The algorithm to find strings requires parallel trans- 
port of the variables at vertex i along a geodesic on CP 2 to 
the vertex j. Then the transported variables are rotated to 
the assigned variables at j, by using an 5(7(2) geodesic trans- 
formation. 



To get a better physical sense for this algorithm, it is 
useful to consider monopole and string formation in the 
simpler symmetry breaking pattern 



SU{2) -> 17(1) 



1. 



(36) 



This example is discussed in Appendix [E] Wc should 
also add that the natural language for our discussion is in 
terms of fiber bundles since what we have in our model is 
an S 3 /Z2 fiber over a CP 2 base manifold. The topology 
of the base manifold, CP 2 , gives rise to monopoles while 
the topology of the fiber, S* 3 /Z2, gives rise to strings that 
may end on monopoles. 



III. RESULTS 

The simulations were done on a cubic lattice of side 
12 i.e. in 24 x 12 3 tetrahedral cells and was repeated 
10 times to gain statistics. The probability of having a 
monopole or antimonopole in a cell is 0.17. If N is the to- 
tal number of string segments, then the relative numbers 
of segments in closed loops, string segments connecting 
like charge monopoles, and string segments connecting 
oppositely charged monopoles, are given by 



Moops 

N 

N±± 
N 

N+- 
N 



0.4%. 
4.2%. 
95.4%. 



(37) 



Sij Rij • 



This shows that roughly 4% of SU (3) monopoles will end 
up in the doubly charged state and survive annihilation 
due to strings. 




5 10 15 20 
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FIG. 6: A cluster of 6 monopoles can form a singlet of SU(3) 
and SU(2), as in ordinary baryons. A bound state of a 
monopole and antimonopole is also possible, as in ordinary 
mesons. The SU(3) charge on a monopole is shown in shades 
of grey (or in color) and the SU(2) charge as a ±. We have 
not shown the U(l) charge. Z3 strings are shown as solid 
lines; Z2 strings as dashed lines. 



IV. DISCUSSION 



FIG. 4: Logarithm of average number density of strings con- 
necting monopoles and antimonopoles versus string length. 
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FIG. 5: Logarithm of average number density of monopole- 
monopole(++) and antimonopole-antimonopole( ) connec- 
tions versus string length. 



The length distribution of +— strings is shown in 
Fig. |4j Denoting the number density of these strings, 
i.e. number of segments divided by the volume (12 3 ), by 

, the least-squares linear fit is 

n+_(0 = (0.46 ± 0.08)e- (0 - 31±ao3)i (38) 

The corresponding distribution of H — h and strings is 

shown in Fig. [5] and the fit is 

n±±(0 = (0.02 ± 0.01)e^ a23±ao7 ) z (39) 



We have studied the formation of monopoles connected 
by strings in an SU(3) model and the results for the 
distribution of monopoles and strings are summarized 
in Sec. IIIII Here we discuss qualitatively how a simi- 
lar analysis in realistic grand unified models would pro- 
ceed. Our experience with SU (3) helps us understand 
and appreciate the difficulties that are likely to be en- 
countered. As an example, consider the minimal grand 
unified model based on a SU (5) symmetry group. The 
symmetry breaking pattern is 

SU(5) -> [SU(3) x SU{2) x U(l)]/Z 3 x Z 2 . (40) 

and, if the non-Abelian magnetic charges are confined, 
the relevant symmetry breakings are 

5(7(3) ->Z 3 , SU(2)^Z 2 . (41) 

The fundamental magnetic monopoles carry SU (3) and 
SU(2) charges in addition to the topological U (1) charge. 
Therefore each monopole will get connected to a Z3 
string and another Z2 string. Then isolated clusters 
of monopoles come in two varieties, similar to known 
baryons and mesons, as shown in Fig. [6] However, a 
likely outcome at formation seems to be that, in addi- 
tion to some isolated baryonic and mesonic clusters, the 
monopolc-string network percolates and we essentially 
obtain one giant structure, such as depicted in Fig. [7J 

It seems hard to explicitly confirm if the network per- 
colates, say by numerical simulation. For example, the 
vacuum manifold at the first stage of symmetry breaking 
is 12 dimensional and it also does not fall into a straight- 
forward category like CP™. Determining the distribution 
of strings is also more complicated since the SU(3) break- 
ing leads to Z3 strings. These problems do not seem 
insurmountable but are hard enough that we have not 
attempted to solve them at the present time. 

If very few baryonic clusters form and instead an infi- 
nite monopole-string network forms, our experience with 
string networks jl5l . liH [ttI . fl^ suggests that the net- 
work energy density scales with time and never comes 
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Using (|A1|) this becomes 



FIG. 7: Drawing of an infinite monopole-string network that 
could result from SU(5) grand unified symmetry breaking. 
The three different shades of circles represent the SU(3) color 
charge and the plus-minus symbols within the circles the 
SU(2) charge. The 17(1) (hypercharge) charge has not been 
shown. The isolated clusters of monopoles have to occur in 
SU{3) and 517(2) singlets. 



to dominate the universe. Processes such as monopole- 
antimonopole annihilation and meson formation could 
dissipate the energy of the network at a rate that is deter- 
mined by the Hubble expansion. However, this scenario 
ignores the process of baryon formation from the net- 
work. Depending on the rate of this process, we could 
still have a monopole over-abundance problem coming 
from the production of baryonic clusters. 
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(SinSkjSml ~ SuSknSmj) / <f S P e Pq 



4-7T 

(z*z J )(d q zlz l + z* k d q zi)(d r z* k z m + z* m d r z n ). (A3) 

The contractions lead to factors such as Z^Z — 1 or else 
similar factors with derivatives, such as Z^dqZ, d q Z^Z, 
or dgZ^drZ. Of the eight terms in (|A3|) . four cancel 
in pairs, and the other four are equal in pairs, yielding 
finally 

Q=— f d 2 S p e pqr (d„Z^ d r Z - d.Z^Z Z^d r Z), (A4) 
2m J 

which, using Eq. (fl~3|) . is precisely Eq. (|14| . 



APPENDIX B: 5(7(3) GEODESIC MATRIX 

Here we will construct the SU (3) matrix Rji such that 

RjiZi — Zj. (Bl) 

There can be many such rotation matrices but we will be 
interested only in the geodesic rotation such that 



exp(iMs), 



(B2) 



where M is a linear combination of ST/ (3) generators and 
s is the geodesic distance between Z t and Zj as given in 
Eq. (EU). 

The procedure we will adopt is to first consider the 
special case when Zi = Z = (0, 0, 1) T . In this case, we 
can find Rjq and the corresponding M. Then we extend 
the result to include the case when Zi is arbitrary. 



1. Zi = Zq case: 



APPENDIX A: TOPOLOGICAL CHARGE 

We wish to show that the two expressions for the topo- 
logical charge, Eqs. (fT4]) and (fl"5]) . arc equivalent. 

The demonstration follows by using the SU(3) identity 

fabcTijT kl T^ n = 2i(S in 8 k jS m i - SuSknSmj)- (Al) 

where T a are ST/(3) generators normalized such that 
tr(T a T b ) = 2S ab , and f a bc are the structure constants 
defined by [T a ,T b ] = 2if abc T c . The above identity is a 
generalization of the better known identity for the SU (2) 
generators a a : 

tabc&ij&kl&mn = ^(^inSkj^ml ~ ^a^kn^mj)- (A2) 

Now, if we choose Z^Z = 1, Eq. HH) can be written 

Q = ^fabcT^T^T,^ J d 2 S p 

e^(z* z .) 9g (^) 9r&n ). 



Now 

Zl = (0,0,1). (B3) 

Let us denote 

Zf = (z u z 2 ,z 3 ), (B4) 

where zi,z%,Z3 are complex numbers and we assume 
z\Zj = 1. 

We wish a matrix M such that 

Zj=exp(iMs)Z . (B5) 

The matrix M is a linear combination of SL7(3) genera- 
tors. However, the generators of the unbroken SU(2) x 
C/(l) sub-group need not be included since they have no 
effect on Zq- So we need only consider M of the form 

/ -iv\ 
M = -iw , (B6) 
\iv* iw* I 
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where v, w are complex numbers. M is normalized using 
tr(M 2 ) = 2 and so \v\ 2 + \w\ 2 = 1. 

We want to find v,w hi terms of z\, Z2,z%. By the 
standard procedure of diagonalizing M or by using the 
formula M 3 = M, one finds 



R j0 



iAIs 



\v\ 2 cos s + \w\ 2 —vw*(l — coss) vs'ms' 
-v*w(l — cos s) \v\ 2 + \w\ 2 cos s w sins I (B7) 
— u*sins —it;* sins coss 

Now we can relate v, w to z%, z^, Z3. We have 

'zi 



I V Sill s 

-*] — I z 2 I = RjoZo = 1 
,z 3 l \ c 



and so, in terms of the parametrization l|7|). 



w = cos </)e IQ , w = sin<Ae I/3 . 



(B8) 



(B9) 



Note that, from Eq. (|27|) . the distance between i?o and 
Zj is s. This shows that the matrix exp(iAfs) is indeed 
the SU(3) transformation (labeled by s) that traces a 
geodesic from Zq to Zj. Note also that because in our 
convention (JT]) the third component of Zj is real, there is 
no need for an extra phase factor here. 

It can also be verified by explicit substitution that one 
may write Rjq in terms of Zq and Zj as 



R 



jo 



{Z Q + Zj){zl + Z]) 



1 + z]z 



' IZjZl (BIO) 



Next, we relax the condition Zi = Zq. 



General Zi case: 



We would like to find Rji such that 



RjiZi 



Zj, 



(Bll) 



where Rji = cxp(iAfs) and s is the geodesic distance 
between arbitrary points Zi and Zj in CP 2 . 

We already know how to construct the matrix Rio as 
in Eq. (|B7[) that rotates from Zq to Z\. Next find the 
point 



(B12) 



where the bar on the subscript j in Zj denotes that the 
point is obtained by rotating Zj. It is important to note 
that the third component of Zj may not be real. In fact, 
since scalar products are unchanged by SU(3) transfor- 
mations, the third components is z\Zj = Z\Zj. 
Next we find Rj such that 



(B13) 



where to use the result in Eq. (|B7[) or (|B10[) requires 
removing the phase factor, i.e., 



Z]Zi 



R ]0 Z - z j~^ 



Z) Zl 



Then it is straightforward to check that 



Z]Zj 

RjiZi — Zj 



\Z] Zi 



where 



Rji - RioRjoR-lo- 



(B14) 



(B15) 



(B16) 



Note that the rotation Rji is the shortest such rotation 
since Rj is the shortest rotation from Zq to Zj . The i?.;o 
transformations in Eq. (|B16|) translate the geodesic path 
from Zq to Zj such that it now goes from Zi to Zj . 

It is also possible to write an explicit formula analogous 
to (|B10[) for Rji. In fact, we have simply to replace Zq 
in that formula by Z\ and Zj by Zj(Z^Zi/\Z^Zi\). 



APPENDIX C: 



CONSTRUCTION OF THE 
MATRIX S. 



The matrix Sji is an SU (2) geodesic rotation that tran- 
forms (Aji, Bji) to (Aj, Bj) at the point Zj on CP 2 (see 
Fig. [3] and Eq. (|3Tj) ). These are the well-known Euler 
rotations e.g. see Section 4.5 in [2lj . 

First we apply the rotations Rj Q to parallel trans- 
port all quantities from Zj to Zq where we know that 
the unbroken SU(2) lies in the 1-2 block of the genera- 
tors. Quantities at Zq will carry a (0) superscript e.g. 



{Af),BfJ) and (Af, Bf 1 ). Then we perform an 517(2) 



1(0) R (o) 



rotation S 



(0) 



1(0) R (0)n 



that rotates (A^' , B^ 



„™.„ . ^ y.,^ ) ZO [Jij ,£>j ). 

There are two such rotations, each of which can be writ- 
ten as 



S 



(0) 



ji 



e in-<r$/2 q 
1 



(CI) 



where cr denotes the three Pauli spin matrices, and ip, 
and 4> are the Euler angles. The angle of rotation, $, is 
given up to a two-fold ambiguity, 



$ 6 + ip 9 

cos — = ± cos cos - , 

2 2 2' 



and 



with 



sin($/2)' 



6- ib . 6 

e-i = cos sin-, 

2 2 

6- ip 9 
e 2 = sm — - — sin-, 

cb + ib 6 
e 3 = sm — 2 — C0S 2' 



(C2) 
(C3) 



(C4) 
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The Eulcr angles </>, ip and 8 can be written in 

(o) h (o) (0) 

ji > u ji ) *"./' 



terms of the vector triads at Zq, (a^° , b,°\ c^) and 



(a^ ,b^- , c^- ) where c = axb: 



cos6> 


= c (0) • 


c (0) 










cos^ 


= „?>■ 


c, 




sin-0 


= (af 


xC) 


(0) 


C0S</> 


- a (0) - 


c 




sin 






■ c (0) 



(C5) 



where £ is a unit vector along the "line of nodes" 



eff x cf 
Ic^xcfl 



Finally, the matrix S^P can be parallel transported 



(C6) 



back to Zj to obtain 



Sji — RjoSjJ Rj . 



(C7) 



The two-fold ambiguity in the rotation corresponds to 
two possible angles of rotation, by $ or by $ — 2ir. We 
choose the rotation that is smaller i.e. |$| < tt. 



APPENDIX D: CONSISTENCY OF MONOPOLE 
AND STRING NUMBERS 

The topology of the symmetry breaking scheme de- 
scribed by Eqs. (|TJ) followed by ([2j> requires that a cell 
with a nonzero monopole number has an odd number of 
strings through its faces, while one with zero charge has 
an even number. Here we demonstrate that the formal- 
ism described above respects this condition. 

For this purpose it is convenient to rotate all the rel- 
evant quantities to the base point Zq. In particular, we 
consider, in place of (|35[) the quantity 



W 



(0) 

{ijk} 



? (0) 



= S{S fe} exp(i*a tofe} ^/3T 8 ), (Dl) 



where 



(D2) 



Clearly, W^j k y must be one of the two central elements 

of SU(2)q, and consequently Sffio € U(l)o since the 
other two factors in (|Dlj) are in that subgroup. 

Now consider the product of the W'°'s from all four 
faces, say 

" "{123} "{142} "{134} ^ {243}" 



The order of the four factors is arbitrary but has been 
chosen for later convenience. This product is evidently 
again one of the two central elements of SU (2)o; which 
one determines whether the number of strings entering 
the cell is even or odd. 

Since T 8 commutes with all the S^ k y, when we substi- 



tute from (|Dlj) into (|D3j) , we can move all the exponential 
factors to the right, and so write as a product 



W {0) = S (0) exp(z^QV3T 8 ), 
where we have used Eq. (f2"Tj) . and 



5 (o) 



c(0) c-(O) c(0) 

Jr-inni U r i I ii l iJ r 



"'{123}^{142}' J {134}' J {243} 



s 



(0) 



(D4) 



(D5) 



Moreover, using Eq. (|32[) . we see that each factor here 
may be written as a product of three factors coming from 
the edges of the triangle, each transported to Zq: 



c-(o) _ n (o) n (o) T Ao) 

^{ijk} ~ U ik U kj U ji ! 



(D6) 



where, for example, 



(D7) 



The key now is to compare the transformations U. 

AO) 



(0) 

ji 



and U>j . By construction, SjiRji transforms $>i, ^u, 
hito $j, &ij, x &2j, whereas SijRij performs the 
inverse transformation. Moreover, the prescription for 
choosing between the two possible transformations is the 
same in each case. These two products are therefore in- 
versses. Thus we learn that 



rr(0)t = rr(0) 

ji ij ' 



(D8) 



Now when we substitute (|D6[) into (|D5[) we find 

C(0) _ 7-r(0)rr(0)rr(0) j M) ,,(0) j M) 

a — u 13 u 32 u 21 .u 12 u 2i u 41 

x rr(0)rr(0) r/ (0) n (0) n (0) TI (0) 
XL/ 14 u 43 u 31 - u 23 u 34 u 42 • 



(D9) 



These factors are six pairs of mutual inverses, although 
since they do not necessarily commute, it is not immedi- 
ately obvious that they cancel. It is clear, however, that 
two pairs cancel at once, leaving us with 

SCO) = ^(0)^(0) ^(0) ^(0)^(0) ^(0)^(0)^0) . (D1Q) 

But now recall that the product of the last three factors 
is S|24 3 i € U(1)q. Consequently, this product commutes 

with all the f/ (0) s, so we may move these three factors 
together to any desired position in the product. Placing 
them after the first two we find 

S(°) = E4 0) t4 0) ■ U^Mf ■ t4 0) • . (Dll) 

But now it is clear that we can cancel these pairs succes- 
sively, so that finally we obtain 



s<°) 



(D12) 
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So this factor may be cancelled from the right side of 
Eq. (|D4[) , which then becomes 

W {0) = cxp(i7rQV3T 8 ), (D13) 

This shows, as required, that the number of strings is 
odd or even according as Q = 1 or 0. 

APPENDIX E: SU(2) MONOPOLES AND 
STRINGS 

Here we discuss monopoles connected by strings in the 
model 

SU(2)->U(1)->1. (El) 

The first symmetry breaking is achieved by giving a VEV 
to an SU(2) adjoint, equivalent to choosing a unit 3- 
vector (call it v). The vacuum manifold is SU(2)/U(1) = 
S 2 . The second symmetry breaking is achieved by giving 
a VEV to a second SU(2) adjoint, call it a, which is 
orthogonal to v. At this stage the vacuum manifold is S 1 . 
Therefore monopoles are formed in the first symmetry 



breaking and these get connected by strings in the second 
symmetry breaking. 

To simulate monopole formation, we assign unit vec- 
tors v, equivalentl y p oints on S 2 , to the points on our 
spatial lattice [H, H3]- A tetrahedral cell gets mapped 
to a tetrahedron in S 2 and some of these mappings will 
be incontractable, implying the existence of a monopole 
within the tetrahedral cell. 

The formation of strings that connect the monopoles 
is more involved but easy to picture, as in Fig. [8] Since a 
is orthogonal to v, we can view it as picking a direction 
on the tangent plane of the S 2 . To determine if there 
is a string passing through a triangular plaquette of the 
spatial lattice, we have to parallel transport a between 
the vertices of the triangle using rotations R and then 
rotate the transported vectors at the vertices using S. 
This is explained in Fig. [8] The scheme for SU(3) is just 
a generalization of the scheme for the SU (2) model. The 
complications are technical in that, instead of the tangent 
plane, our "vectors" at every vertex lie on an S 3 /Zq fiber 
and the geodesies and rotations are harder to determine 
in practice. 
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FIG. 8: To determine if a string passes through a spatial 
triangular plaquette, we first take the corresponding triangle 
on S 2 , labelled {ijk}, and then determine if the vector in the 
tangent plane rotates by 2ir in circumnavigating the spherical 
triangle. To do this, we first parallel transport the vector from 
i to j along a geodesic, described here as a rotation, R. Then 
we find the rotation S within the tangent plane that takes the 
transported vector into the vector at the vertex j. In each case 
we choose the minimal-angle rotation. Then we do the same 
thing for the remaining sides. Since we end up with the same 
vector at i that we started with, the combined transformation 
is either the identity or a 2n rotation. 



